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ABSTRACT 


Cross-shore  exchange  on  the  inner  shelf  has  important  impacts  on  the 
ecosystem,  transporting  heat,  nutrients,  pollutants  and  phytoplankton  between 
the  midshelf  and  surf  zone.  The  effects  of  a  strong  (cross-shore  wind  stress,  Tsx 
>0.05Pa)  diurnal  (7-25  hrs)  sea  breeze  on  cross-shore  exchange  at  Marina, 
Monterey  Bay,  California  is  investigated  using  two  years  of  continuous  winds, 
waves,  and  ocean  velocities.  Surface  wind  stress  has  spectral  peaks  at  1,  2,  and 
3  cpd  and  the  diurnal  wind  variability  is  greater  than  50%.  Similar  spectral 
energetic  peaks  also  occur  with  waves  and  currents.  During  sea  breeze 
relaxation  (-0.05Pa  <rsx<0.05Pa),  a  background  wave-driven  inner-shelf 
undertow  profile  exists,  which  is  equal  and  opposite  to  Lagrangian  Stokes  drift, 
resulting  in  a  net  zero  Lagrangian  transport  at  depth.  In  the  presence  of  a  sea 
breeze  {Tsx  >0.05Pa),  a  uniform  offshore  profile  develops  that  is  different  from  the 
background  undertow  profile  allowing  cross-shore  Lagrangian  transport  to 
develop,  while  including  Lagrangian  Stokes  drift.  The  seasonality  of  waves  and 
winds  modify  the  diurnal  sea  breeze  impact.  Therefore,  material  is  hypothesized 
to  incrementally  move  onshore  near  the  surface  and  offshore  near  the  sea  bed 
only  during  sea  breeze  events. 
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I.  INTRODUCTION 


Cross-shore  transport  plays  a  significant  role  in  the  allocation  and 
redistribution  of  phytoplankton,  nutrients,  sediment  and  pollutants  across  the 
continental  shelf  and  the  surf  zone.  Continental  shelf  ecosystems  are  some  of 
the  most  dynamic  on  earth  (Falkowski  et  al.  1998).  Near  the  continental  shelf 
break,  upwelling  forces  nutrient-rich  water  to  the  surface  causing  high  biological 
productivity,  which  is  then  transported  across  the  continental  shelf  (Pilskain  et  al. 
1996).  Therefore,  it  is  important  to  understand  the  mechanisms  responsible  for 
cross-shore  exchange.  There  are  three  primary  regions  located  shoreward  of  the 
outer  continental  shelf,  which  are  referred  to  as  the  midshelf,  inner-shelf  and  surf 
zone,  and  are  dominated  by  different  dynamics  (Lentz  1999). 

For  the  midshelf,  cross-shore  transport  is  the  result  of  steady  alongshore 
wind  stresses  acting  together  with  the  Coriolis  force  producing  a  transport  of 
water  at  a  right  angle  to  the  predominant  wind  direction  in  the  Northern 
Flemisphere,  known  as  Ekman  transport  (Ekman  1905,  Sverdrup  1938,  Csanady 
1978,  Lentz  1992,  Dever  1997a).  On  the  midshelf,  during  upwelling  or 
downwelling  events,  the  midwater  pycnocline  intersects  the  surface  or  bottom 
boundary,  forming  a  front  that  moves  offshore.  During  both  upwelling  and 
downwelling,  full  Ekman  transport  develops  in  the  midshelf  region  because  the 
strong  stratification  of  the  pycnocline  acts  to  insulate  the  surface  and  bottom 
layers  from  each  other  (Austin  and  Lentz  2002).  The  region  seaward  of  this 
dynamic  front  is  the  midshelf  and  the  region  shoreward  of  the  front  is  the  inner- 
shelf  (Lentz  1994,  Austin  and  Lentz  2002). 

On  the  inner  shelf,  as  the  water  depth  decreases,  the  alongshore  surface 
stress  becomes  increasingly  balanced  by  the  bottom  stress  instead  of  the 
Coriolis  force,  reducing  Ekman  surface  boundary  layer  transport  (Dever  et  al. 
2006).  Numerical  model  studies  (Austin  and  Lentz  2002,  Tilburg  2003)  and 
observations  on  the  inner-shelf  regions  of  North  Carolina  (Lentz  1999,  2001), 
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Massachusetts  (Fewings  et  al.  2008),  Oregon  (Kirincich  et  al.  2005)  and 
California  (Cudaback  et  al.  2005)  found  that  alongshore  winds  are  not  a  sufficient 
mechanism  in  driving  cross-shore  circulation  on  the  inner-shelf.  Tilburg  (2003) 
numerically  found  that  cross-shore  winds  induced  a  two-layer  circulation 
consisting  of  onshore  transport  near  the  sea  surface  and  an  equal  and  opposite 
offshore  transport  below  that  allowed  for  cross-shore  exchange. 


The  surf  zone  is  defined  as  the  region  between  depth-limited  breaking 
(Thornton  and  Guza,  1983)  and  the  shoreline.  The  mechanisms  for  cross-shore 
transport  in  the  surf  zone  are  wave-driven  and  consist  of  Stokes  drift,  undertow, 
and  rip  currents.  Stokes  drift  is  the  time-averaged  second  order  velocity  of  a 
fluid  particle  under  a  wave.  There  is  an  incomplete  closure  of  the  parcel  path 
after  each  wave  period  resulting  in  a  net  drift  in  the  direction  of  wave  propagation 
(Stokes  1847,  Madsen  1978).  The  associated  mass  transport  occurs  between 
the  wave  trough  and  crest  in  the  Eulerian  frame  and  is  vertically  distributed  below 
mean  sea  level  (MSL)  in  the  Lagrangian  frame.  The  theoretical  cross-shore 
Stokes  transport,  is. 


a 


stokes 


SK  n 

- ^cos6'„ 

16c 


(1) 


where  g  is  the  gravitational  acceleration,  7/^,^  is  the  significant  wave  height,  cis 


the  phase  speed  of  the  waves,  and  0^  is  the  wave  direction  relative  to  shore 
normal  (Stokes  1847,  Longuet-Higgins  1953).  The  onshore  transport  near  the 
surface  is  balanced  by  an  equal  transport  in  the  opposite  direction  at  depth,  the 
undertow  (Ursell  1950,  Haines  and  Sallenger  1994,  Reniers  et  al.  2004).  In  the 
surf  zone,  the  compensating  return  2D  vertical  profile  is  parabolic  with  a 
maximum  offshore  flow  at  mid-depth  and  close  to  wave  breaking. 


Undertow  affects  the  2D  vertical  profile,  while  rip  channels  affect  the  2D 
horizontal  circulation  (assuming  the  vertical  is  depth-uniform)  allowing  for 
transport  of  material  across  the  surf  zone.  Alongshore  bathymetric  variations 
create  alongshore  variations  in  depth-limited  wave  breaking  that  induce 
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alongshore  gradients  in  pressure  and  momentum  driving  water  from  the  shore- 
connected  shoals  toward  rip  channels  (Bowen  1969,  Dalrymple  1978).  Onshore 
flows  over  the  shore-connected  shoals  or  bars  transition  to  alongshore  flows 
(feeder  currents)  near  the  shoreline  that  converge  in  the  rip  channel  and  extend 
seaward  across  the  surf  zone,  and  beyond  (Shepard  et  al.  1941,  Inman  and 
Brush  1973).  Inman  and  Brush  (1973)  referred  to  rip  currents  as  “freeways”  that 
transport  material  outside  the  surf  zone  to  the  inner-shelf. 

A.  MECHANISMS  FOR  CROSS-SHORE  TRANSPORT  ON  THE  INNER- 


SHELF 


Hasselmann  (1970)  suggested  that  the  Coriolis  force  acting  on  the  surface 
wave  flow  induces  a  small  along-crest  wave  velocity  that  is  in  phase  with  the 
vertical  wave  velocity  resulting  in  a  “wave  stress”.  Xu  and  Bowen  (1994) 
theoretically  determined  that  the  “Hasselmann”  wave  stress  is  balanced  by  the 
Coriolis  force  in  the  alongshore  momentum  balance,  resulting  in  an  offshore- 
directed  Eulerian  flow,  which  is  equal  and  opposite  to  the  Lagrangian  Stokes 
drift,  resulting  in  a  zero  net  transport  over  the  vertical  as  suggested  by  Ursell 
(1950).  This  would  explain  why  the  ocean  does  not  pile  up  with  water  at  the 
shoreline  owing  to  the  shoreward  transport  of  Stokes  drift.  The  presence  of  this 
compensating  offshore-directed  Eulerian  flow  has  been  observed  outside  of  the 
surf  zone  in  the  field  (Reniers  et  al.  2004,  Smith  2006,  Lentz  et  al.  2008)  and  in 
laboratory  studies  (Putrevu  and  Svendsen  1993,  Ting  and  Kirby  1994,  Cox  and 
Kobayashi  1997,  Monosmith  et  al.  2007). 

Lentz  et  al.  (2008)  compared  the  wave-driven  velocity  profiles  from 
Martha’s  Vineyard  Coastal  Observatory,  MVCO,  to  the  theoretical  work  by  Xu 
and  Bowen  (1994).  The  model  by  Lentz  et  al.  (2008)  assumes  steady-state, 
linear  dynamics  with  alongshore  homogenous  variations,  and  constant  density. 
The  model  is  based  on  the  continuity  and  momentum  balances. 


(2) 
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-fv  =  fv^  -  gri^  -  F'*'^  +  I +(Au^)z  > 
fu  =  -/u^f  +  T^^^'/Po  +  (^Vz)z  - 


(3) 

(4) 


where  Q^tokes  is  the  onshore  Stokes  transport,  fv  and  fu  are  the  Coriolis  forces, 
fv^^  and  fu^j.  are  the  Hasselmann  wave  stresses,  g  is  gravitational  acceleration, 
r}^  is  the  cross-shelf  pressure  gradient,  F'*’®  is  momentum  flux  divergence, 

Ip^  and  t^'^^VPo  are  the  near  bottom  wave  stresses,  and  are 

the  turbulent  eddy  viscosities  (see  Lentz  et  al.  (2008)  for  a  detailed  description  of 
the  model).  For  small  eddy  viscosities  (indicative  of  motions  outside  the  surf 
zone)  the  offshore  velocity  profiles  had  a  curvature  with  a  maximum  offshore  flow 
near  the  surface,  that  decreases  with  decreasing  depth.  These  modeled  profiles 
favorably  matched  the  field  observations  at  MVCO  outside  the  surf  zone.  In  the 
absence  of  wind,  Lentz  et  al.  (2008)  concluded  that  the  time-averaged  flow  in  the 
inner  shelf  is  primarily  associated  with  undertow  that  is  driven  by  surface  gravity 
waves  and  influenced  by  the  Hasselmann  wave  stress. 

Prior  to  the  work  of  Tilburg  (2003),  cross-shelf  winds  were  considered 
ineffective  in  forcing  cross-shelf  transport  due  to  the  cross-shelf  components  of 
the  surface  and  bottom  boundary  stress  being  an  order  of  magnitude  smaller 
than  the  Coriolis  force  of  the  alongshelf  flow  (Csanady  1978;  Allen  1980). 
Fewings  et  al.  (2008)  observed  that  the  cross-shore  wind  forcing  was  important 
for  cross-shelf  exchange,  as  it  modifies  the  wave  driven  undertow  profile  found 
by  Lentz  et  al.  (2008),  inducing  a  two-layer  flow  circulation.  The  surface  flow 
below  the  wave  trough  was  in  the  direction  of  wind  forcing,  consistent  with 
previous  estimates  of  surface  wind-induced  drift  (Wu,  1983,  Ogasawara  and 
Yasuda  2004),  and  a  compensating  return  flow  in  the  bottom  portion  of  the  profile 
was  in  the  opposite  direction  of  the  wind  forcing.  The  presence  of  a  cross-shore 
wind  results  in  a  non-zero  net  transport  at  depth.  Fewings  et  al.  (2008)  found 
that  as  background  wave  forcing  increased,  the  two-layer  circulation  resulting 
from  wind  forcing  was  subjugated  by  the  wave-driven  undertow,  which  is  always 
directed  offshore.  When  wave  and  wind  forcing  were  in  the  same  direction,  the 
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cross-shelf  velocity  profile  is  practically  vertically  uniform  owing  to  the  vertical 
shears  being  similar  in  magnitude  but  opposite  in  direction.  When  wind  and 
wave  forcing  were  in  opposing  directions,  the  cross-shelf  velocity  profile  is 
strongly  sheared  owing  to  the  vertical  shears  being  large  and  in  the  same 
direction. 

The  Lagrangian  velocity,  ul,  is  defined  as, 

Ml  =«£+«.  (5) 

where  Ue  is  Eulerian  velocity  and  us  is  Stokes  drift  velocity.  For  the  wave-driven 
undertow  on  the  inner-shelf,  Ue  and  us  are  equal  and  opposite  at  depth  resulting 
in  zero  ul  at  depth  (Xu  and  Bowen,  1994;  Lentz  et  al.  2008).  This  results  in  zero 
net  transport.  However,  when  the  undertow  profile  is  modified  due  to  winds 
(Fewings  et  al.  2008),  the  transport  is  onshore  (offshore)  near  the  surface  and 
offshore  (onshore)  near  the  sea  bed  with  onshore  (offshore)  winds,  with  the 
depth-integrated  transport  equal  to  zero. 

Lentz  et  al.  (2008)  confirmed  that  the  forcing  mechanisms  for  cross-shelf 
exchange  on  the  inner  shelf  are  fundamentally  different  from  those  at  midshelf. 
Fewings  et  al.  (2008)  evaluated  the  relative  importance  of  cross-  and  alongshore 
wind  and  wave  forcing  as  a  function  of  water  depth  and  suggested  that  cross¬ 
shore  wind  stress  is  the  primary  forcing  mechanism  for  cross-shelf  exchange  for 
water  depths  less  than  30  m. 

B.  SEA  BREEZE 

The  diurnal  sea  breeze  is  the  result  of  a  difference  in  air  pressure  between 
the  land  and  ocean  caused  by  differential  heating,  which  occurs  for 
approximately  two-thirds  of  the  earth’s  coasts  (Defant  1951,  Atkinson  1981, 
Simpson  1994).  Owing  to  a  difference  in  thermal  heat  capacity,  the  land  heats 
and  cools  more  rapidly  than  the  ocean,  generating  an  intense  onshore  sea 
breeze  in  the  late  afternoon  and  a  weaker  offshore  land  breeze  during  the  night 
(Defant  1951,  Atkinson  1981,  Simpson  1994). 
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Woodson  et  al.  (2007)  observed  that  the  diurnal  sea  breeze  in  northern 
Monterey  Bay  was  responsible  for  local  upwelling  resulting  in  offshore  Ekman 
transport  during  the  day.  The  sea  breeze  at  this  site  is  predominantly  parallel  to 
the  shoreline  (Figure  1).  Note  that  our  study  site  is  located  in  the  center  of 
Monterey  Bay  and  focuses  on  similar  events  (Figure  1),  where  the  winds  are 
predominantly  perpendicular  to  the  shoreline.  Woodson  et  al.  (2007)  suggests 
that  this  local  diurnal  upwelling  may  supply  nutrients  to  nearshore  kelp  beds  and 
transport  larvae  to  nearshore  habitats. 

The  effects  of  a  diurnal  sea  breeze  and  waves  on  cross-shelf  exchange  on 
the  inner  shelf  are  examined  at  Marina,  Monterey  Bay,  CA.  It  is  hypothesized 
that  the  sea  breeze  diurnally  modifies  the  vertical  current  profile  allowing  for 
cross-shore  exchange  to  develop.  Two  years  of  continuous  ocean  velocity, 
wave,  and  wind  data  are  analyzed  and  described  below  for  testing  the  proposed 
hypothesis. 
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II.  OBSERVATIONS 


The  subaqueous  ocean  velocity  profiles,  with  a  bin  size  of  0.5  m,  and  co¬ 
located  pressure  were  measured  continuously  at  1  Hz  by  a  bottom-mounted 
1200-kHz  broadband  Acoustic  Doppler  Current  Profiler  (ADCP)  deployed  in 
approximately  13  m  water  depth  at  36°42.0’N,  121°48.9’W  (Figure  1)  from  01 
January  2006  to  31  December  2008.  Surface  winds  and  air  temperature  were 
measured  every  2  minutes  by  a  wind  anemometer  mounted  on  a  10  m  tripod  on 
a  local  sand  dune  that  was  51  m  above  sea  level  and  approximately  2.5  km  east 
of  the  ADCP  site  at  Marina  Airport  (36°40.9’N,  121°46.0’W)  (Figure  1).  Daily 
records  with  greater  than  30  percent  erroneous  or  missing  data  were  not 
considered  in  the  analyses.  The  remaining  data  were  quality-controlled  by 
applying  a  3-standard  deviation  filter  to  remove  outliers  (Emory  and  Thomson, 
2001). 

Since  it  is  difficult  to  recognize  the  sea  surface  with  ADCP  velocity 
observations,  sea  surface  elevations  were  estimated  by  transforming  the  co¬ 
located  ADCP  pressure  measurements  applying  linear  wave  theory.  Previous 
studies  have  not  used  the  ADCP  velocity  measurements  above  the  wave  trough 
due  to  acoustic  surface  reflections  (Fewings  et  al.  2008,  Lentz  et  al.  2008)  and 
difficulties  identifying  the  sea  surface.  Instead,  wave-driven  surface  transport 
(Stokes  drift)  is  estimated  using  linear  wave  theory.  The  variability  of  surface 
transport  by  waves  and  winds  is  an  important  aspect  of  this  research,  so  ADCP 
velocities  above  the  trough  are  included.  ADCP  velocities  above  the  calculated 
sea  surface  elevation  were  set  to  zero.  The  1  Hz  velocity  profile  was  depth- 
normalized  by  using  a  low-pass  filtered  (/<0.04Hz)  sea  surface  elevation,  which 
includes  infragravity  waves  (Herbers  et  al.  1994)  and  tides.  This  allows  the  high- 
frequency  waves  (/>0.04  Hz)  associated  with  sea  and  swell  waves  to  fluctuate 
around  z/h=1,  referred  to  as  mean  sea  level  (MSL).  The  depth  normalization  also 
removes  the  tidal  modulation  allowing  the  time  records  greater  than  3  hours  to  be 

compared.  The  normalized  scale  was  set  from  z/h  =  0  at  the  sea  bed  to  z/h=  ^  .5 
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at  the  highest  wave  crest  with  0.05  vertical  increments.  For  the  purposes  of 
calculating  transport,  the  flow  was  considered  vertically  uniform  between  the  first 
good  bottom  bin  {z/h  =  0.2)  and  the  sea  bed. 

Cross-shore  velocity  magnitudes  are  easily  affected  by  slight  rotational 
errors  in  the  local  coordinate  frame  due  to  the  tidal  currents  being  predominantly 
alongshore  (Rosenfeld  et  al.  2008).  A  Principal  Component  Analysis  (PCA)  was 
used  to  determine  the  local  coordinate  frame  orientation  (Figure  2),  where  the 
alongshore  velocity  component  is  defined  as  the  major  axis  of  the  subtidal  depth- 
averaged  velocity,  Uda  ,  computed  from  zero  depth-crossing  to  the  sea  bed 
(Davis,  1986).  For  quality  control  purposes,  the  PCA  was  estimated  for  all 
normalize  depths  to  prevent  biases  associated  with  surface  gravity  waves  and 
bottom  drag  owing  to  to  bed  stress.  It  was  found  that  between  z/h  =  0.25  and  z/h 
=  0.75,  the  PCA’s  were  within  1°.  Below  z/h  =  0.25  and  above  z/h  =  0.75  the 
PCA  rotational  angles  were  inconsistent.  The  final  PCA  used  a  Uda  between  the 
depths  of  z/h  =  0.25  and  z/h  =  0.75  indicating  that  the  major  principal  axis  was 
oriented  10°  clockwise  from  true  north,  consistent  with  local  isobaths  from  a 
recent  nautical  chart  (  Monterey  Bay  18685  33'^'^  edition  September  2005)  (Figure 
1).  A  right-handed  coordinate  system  is  used,  where  x  is  positive  cross-shore 
eastward,  y  is  positive  alongshore  northward,  and  z  is  positive  upward. 

At  this  location,  the  winds  are  predominantly  cross-shore  making  it  an 
ideal  location  for  studying  cross-shore  sea  breezes  effects.  Surface  wind 
stresses  {Ts)  were  calculated  using  the  model  from  STRESSTC  (Smith  1988) 
using  wind  speed  and  air  temperature  at  a  height  of  10  m  above  ground.  Cross¬ 
shore  surface  wind  stress  (Tsx)  and  along-shelf  surface  wind  stress  (fsy  )  were 
rotated  to  the  same  local  coordinate  system  defined  by  Uda  The  local  wind  stress 
rose  histogram  indicates  that  90%  of  the  stresses  resided  between  +/-  20 
degrees  of  shore-normal  with  a  mean  of  0.074  Pa  and  variance  of  0.024  Pa^ 
(Figure  3).  The  mean  wind  direction  is  +  6  degrees  relative  to  shore-normal.  The 
diurnal  Tsx  '\s  skewed  (skewness=1 .61 )  onshore,  fsxhas  intense  onshore  stress 
that  last  for  a  short  duration  followed  by  a  longer  minimal  offshore  stress.  Unlike 

8 


the  MVCO  field  site  (Fewings  et  al.  2008),  Monterey  Bay  has  a  strong  and 
consistent  year  round  diurnal  sea  breeze  (Figure  4). 


The  inner  shelf  was  previously  defined  as  the  region  inshore  of  the 
upwelling  and  downwelling  front  where  the  surface  and  bottom  Ekman  layers 
interact  (Lentz  1994,  Austin  And  Lentz  2002,  Tilburg  2003).  Tilburg  (2003)  uses 
the  formula  presented  in  Weatherly  and  Martin  (1978)  to  determine  whether  the 
site  is  located  in  the  midshelf  or  the  inner  shelf. 


S  = 


1 


where  is  a  fitting  parameter  and  is  estimated  by  Tilburg  to  be  2.3, 


(6) 

f  is  the 


Coriolis  parameter  and  N  is  the  buoyancy  frequency.  When  5  is  greater  than 


the  local  water  depth  (/?),  the  region  is  considered  to  be  the  inner  shelf.  The 
average  value  for  N  in  Monterey  Bay  is  3.5  x  10'^  s'^  with  an  /"  =  8.7  x  10'^  s'^ 

(Petruncio  et  al.  1998,  Rosenfeld  et  al.  1994).  At  the  measurement  site  in 
Monterey  Bay  ^  0.005  Pa  99%  of  the  time  resulting  in  5  greater  than  h 
suggesting  this  site  is  always  in  the  inner  shelf.  Though  stratification  may  be 
important  at  different  times  of  the  year,  it  is  assumed  small  for  this  analysis. 
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III.  RESULTS 


A.  SPECTRAL  ANALYSIS 

Record  composed  of  20-minute  block  averaged  values  over  730  days. 

Spectra  were  calculated  using  a  Hanning  window  with  50%  overlap  resulting  in 

102  degrees  of  freedom  and  frequency  resolution  of  0.003  cycles  per  hour  (cph). 

Spectra  are  examined  in  the  diurnal  frequency  band  (0.03-0.13  cph).  Energetic 

spectral  peaks  for  wind  stresses  and  waves  occur  at  diurnal,  1  cycle  per  day 

(cpd)  (/  =  0.042  cph),  2  cpd,  semi-diurnal,  (/  =  0.083  cph)  and  3  cpd  (  /  =  0.125 

cph)  (Figure  5).  The  energy  at  1  cpd  peak  is  the  greatest  and  is  three  times 

larger  than  at  2  cpd  and  ten  times  larger  than  at  3  cpd,  where  the  2  and  3  cpd 

peaks  are  the  wind  harmonics.  Velocity  spectra  were  estimated  at  each 

elevation  and  vertically  stacked.  The  2-D  cross-shore  velocity  spectra  indicates 

a  significant  diurnal  peak  throughout  most  of  the  column  and  less  energetic 

peaks  near  the  surface  at  2  and  3  cpd  (Figure  5).  The  2-D  alongshore  velocity 

spectra  are  an  order  of  magnitude  greater  than  in  the  cross-shore  and  have 

broader  energy  at  1  and  2  cpd  that  occur  throughout  the  column.  The  energetic 

spectral  peaks  are  associated  with  tidal  constituents  (i.e.,  SI,  K1,  M2,  S2,  and 

S4)  and  the  diurnal  wind  harmonics  (Militello  and  Kraus  2001).  There  is  a 

significant  coherence  at  the  95%  level  between  the  cross-shore  wind  and  the 

cross-shore  velocity  at  1,  2,  and  3  cpd  (not  shown).  There  is  no  significant 

coherence  between  the  cross-shore  winds  and  the  alongshore  velocities  for  the 

same  frequencies  (not  shown).  In  order  to  resolve  the  K1  and  1  cpd  spectral 

energy  peaks,  a  maximum  frequency  resolution,  Af,  must  be  defined, 

1  1 

The  window  length  is  1/Af.  For  this  example,  a  minimum  of  1367  days  (3.75 
years)  is  required,  while  only  providing  2  DOF.  The  observations  described 
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herein  for  2  years  shorter  than  the  3.75-year  requirement.  Therefore,  spectral 
analysis  is  unable  to  resolve  these  peaks. 

B.  DIURNAL  VARIABILITY 

Diurnal  variability  is  defined  as  the  variance  of  the  high-pass  (<25hrs) 
filtered  signal  divided  by  the  variance  of  the  total  signal.  The  cross-shore  sea 
breeze  represents  over  56%  of  the  total  surface  wind  variability  (Table  1).  The 
diurnal  variability  of  waves  is  about  20%,  which  is  influenced  by  the  diurnal  sea 
breeze.  The  diurnal  variability  of  the  surface  velocities  {z/h  =  1)  is  24%  in  the 
cross-shore  and  28%  in  the  alongshore,  whereas  below  the  trough,  the  Uda 
diurnal  variability  is  62%  in  the  cross-shore  and  32%  in  the  alongshore. 

Since  many  of  the  significant  tidal  constituents  occur  within  the  diurnal 
frequency  band,  tidal  current  variability  is  examined  (Table  2).  The  tidal  current 
variability  is  calculated  as  the  variance  of  the  predicted  diurnal  tidal  currents 
found  using  the  T_TIDE  program  (Pawlowicz  et  al.  2002),  divided  by  the  variance 
of  the  total  current  at  each  z/h.  Previous  studies  suggest  that  a  record  of  366  or 
more  days  is  needed  to  resolve  S1  from  K1  using  T_TIDE  (Rosenfeld  et  al.  2008, 
Petruncio  et  al.  1998).  Variability  of  the  calculated  depth-averaged  tidal  velocity 
in  the  cross-  and  alongshore  are  9%  and  20%  (Table  1).  However,  the  2  and  3 
cpd  sea  breeze  harmonics  coincide  with  the  S2  and  S4  tidal  constituents  and 
cannot  be  separated  using  T_TIDE.  Militello  &  Kraus  (2001)  found  it  difficult  to 
differentiate  diurnal  tidal  constituents  from  wind  harmonics.  It  is  assumed  that  the 
mean  sea  surface  elevation  is  not  influenced  by  the  sea  breeze,  whereas  the 
velocities  are  influenced.  T_TIDE  applied  to  the  averaged  sea  surface  elevation 
predicts  99.9%  of  the  observed  signal.  The  predicted  SI,  S2,  and  S4  combined 
tidal  elevation  represents  0.5%  of  the  total  prediction.  Currently,  we  are  unaware 
of  any  Fourier  method  that  exists  to  distinguish  motions  at  the  same  frequency 
with  different  phases  (Appendix  A).  It  is  assumed  that  the  velocities  associated 
with  S2  and  S4  are  weak,  owing  to  the  less  than  0.5%  surface  contribution,  and 
predominantly  in  the  alongshore  (Rosenfeld  et  al.  2008). 
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Owing  to  the  small  cross-shore  velocities,  a  typical  day  is  divided  into  24 
one-hour  increments,  where  a  mean  velocity  for  each  hour  in  the  day  is  averaged 
over  the  2  years  (defined  as  the  24-hr  time-averaged  method).  This  approach 
increases  our  confidence  in  the  mean  hourly  velocities  (Bendat  and  Piersol  2000) 
that  varies  over  the  day,  which  is  hypothesized  to  be  related  to  the  diurnal  sea 
breeze.  This  approach  retains  the  background  flow  field  at  this  site.  The 
confidence  error  (±  1 .65  a  /  Vn)  on  the  mean  at  the  90%  significance  level  for  the 
subaqueous  (below  the  wave  trough)  where  a  is  the  standard  deviation  of  the 
velocity  are  0.001  m/s  and  O.OOSm/s  at  the  surface.  The  24-hr  time-averaged 
method  applied  to  the  predicted  tidal  cross-shore  tidal  velocities  for  the 
subaqueous  and  surface  regions  are  0.001  and  0.003  m/s.  Therefore,  the  total 
expected  error  for  the  subaqueous  and  surface  regions  are  0.002  and  0.006m/s, 
where  the  signals  are  an  order  of  magnitude  greater  than  these  errors. 
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IV.  DISCUSSION 


A.  SUBTIDAL  (>25HRS)  MEAN  CURRENT  PROFILES 

Consistent  with  Fewings  et  al.  (2008),  the  20-minute  time-averaged 
velocities  were  de-tided  using  T_TIDE  (Pawlowicz  et  al.  2002)  and  low-pass 
filtered  (>  25  hour)  to  evaluate  the  subtidal  (>25  hrs)  profile.  The  2  years  of 
cross-shore  vertical  velocity  profiles  were  sub-divided  into  four  sub-categories: 
when  both  Tsx  and  Hsig  are  small,  when  either  Tsx  or  Hsig  are  large  and  the  other 
small,  and  when  both  forcing  mechanisms  are  large.  A  one  standard  deviation 
criteria  was  used  to  define  when  Tsx  and  Hsig  were  large  or  small.  Tsx  was 
considered  small  when  less  than  0.05  Pa  and  large  when  greater  than  0.1  Pa. 
Hsig  was  considered  small  when  less  than  0.7  m  and  large  when  greater  than  1.7 
m. 

When  Tsx  and  Hsig  were  small,  the  mean  cross-shore  profile  had  a  small 
offshore  velocity  of  1.5  cm/s  just  below  the  trough  that  diminished  to  zero  half 
way  down  the  column  (Figure  6a).  This  mean  background  cross-shelf  velocity 
profile  was  subtracted  from  the  other  three  mean  profiles  shown  in  Figure  6  b-d. 
During  large  Tsx  (when  Hsig  is  weak),  there  is  an  increase  in  the  onshore  surface 
flow  of  2  cm/s  due  to  wind-induced  drift  (Figure  6b).  The  subaqueous  flow  has  a 
two-layer  structure  with  onshore  flow  of  2  cm/s  in  the  upper  portion  of  the  water 
column  and  offshore  flow  of  1  cm/s  in  the  lower  portion  of  the  column.  These 
cross-shore  results  are  consistent  with  observations  at  MVCO  (Fewings  et  al. 
2008),  Santa  Barbara  Channel  (Cudaback  et  al.  2005)  and  model  results  (Tilburg 
2003).  During  large  Hsig  (when  the  Tsx  is  weak),  the  surface  flow  increased 
significantly  to  1 1  cm/s,  likely  due  to  Stokes  drift  (Figure  6c).  The  subaqueous 
cross-shore  flow  was  offshore  throughout  the  water  column  and  intensified  near 
the  surface  to  1.5  cm/s.  These  observations  are  similar  to  the  inner-shelf 
undertow  profile  by  Lentz  et  al.  (2008).  When  both  Tsx  and  Hsig  are  large,  the 
profile  became  more  vertically  uniform,  consistent  with  Fewings  et  al.  (2008) 
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results.  The  surface  flow  was  onshore  at  16  cm/s  and  below  the  trough  the 
offshore  flow  was  2  cm/s.  In  summary,  the  subtidal  cross-shore  profiles  at 
Monterey  are  similar  to  the  results  by  Lentz  et  al.  (2008)  and  Fewings  et  al. 
(2008).  It  is  important  to  note  that  Tsx  >0.05  Pa  is  the  threshold  considered  to 
modify  the  wave-driven  undertow  profile. 

B.  DIURNAL  (7-25HRS)  VARIABILITY 

The  thermally  induced  sea  breeze  begins  to  develop  at  1630  GMT  in  the 
presence  of  near  zero  background  Tsx  (Figure  7a).  The  Tsx  peaks  at  2200  GMT 
with  a  magnitude  0.08  Pa,  and  exceeds  the  threshold  of  0.05  Pa  suggesting  that 
it  will  modify  the  undertow  profile.  Since  the  sea  breeze  initiates  near  the 
shoreline,  there  is  no  immediate  diurnal  impact  on  Hsig  due  to  a  minimal  fetch 
(Massel  1996).  As  the  sea  breeze  develops  {Tsx  >  0.05  Pa),  its  extends  both 
landward  and  seaward  increasing  the  fetch.  As  the  fetch  increases  (and 
duration),  the  sea  breeze  begins  to  influence  the  waves  at  2300  GMT  and 
peaking  at  0200  GMT  with  a  Hsig  of  1 .4  m  (Figure  7b).  The  sea  breeze  enhances 
Hsig  by  an  average  of  25  cm  or  16  %  of  the  mean  Hsig.  Surface  flow  is  onshore 
throughout  the  day  due  to  a  combination  of  Stokes  drift  and  wind-induced  drift 
(Figure  7c).  As  the  sun  sets,  the  thermal  gradient  dissipates  due  to  the  cooling 
of  the  coastal  land  (Santa  et  al.  1993,  Santa  1995).  The  “typical  day”  in 
Monterey  indicates  no  offshore  flow  from  a  land  breeze  consistent  with  Santa  et 
al.  (1993)  and  Santa  (1995).  This  time  period  of  minimal  diurnal  forcing  (-0.05  Pa 
<  Tsx  <  0.05  Pa)  is  called  the  relaxation  period,  and  is  below  the  threshold  values 
for  wind  to  modify  the  undertow  profile.  Selow  the  trough  during  the  relaxation 
period,  though  winds  and  waves  are  present,  it  is  considered  to  be  wave- 
dominated  with  a  maximum  offshore  of  3  cm/s  near  the  surface,  which  is 
consistent  with  the  Lentz  et  al.  (2008).  These  results  suggest  that  the  undertow 
is  always  present  at  this  site  owing  to  the  background  presence  of  waves.  Once 
the  sea  breeze  develops  {Tsx  >  0.05  Pa),  the  below  trough  profile  is  near  vertically 
uniform  offshore,  similar  to  the  results  for  combined  cross-shore  wind  and  wave 
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forcing  found  in  the  subtidal  results  (Section  4.1)  and  by  Fewings  et  al.  (2008). 
The  background  undertow  profile  is  modified  during  the  sea  breeze  allowing  for 
the  onshore  transport  near  the  surface  and  offshore  transport  near  the  sea  bed, 
which  is  not  reversed  during  relaxation. 

C.  SEASONAL  MODULATION  OF  THE  DIURNAL  VARIABILITY 

The  sea  breeze  persists  year  round  at  Monterey  Bay  and  only  changes 
seasonally  in  intensity  (Figure  4).  The  24  hr  time-averaged  method  was 
evaluated  seasonally  using  the  four  northern  hemisphere  meteorological 
seasons,  Spring  (March  1®*  through  May  31  ),  Summer  (June  1®‘  through  August 
31  ),  Fall  (September  1®‘  through  November  30*^  )  and  Winter  (December  1®* 

through  February  28**^  )  (Trenberth,  1983).  For  seasonal  evaluation,  the 
subaqueous  current  data  were  high-pass  filtered  (<25hr)  in  order  to  enhance  the 
diurnal  variability. 

The  Tsx  and  Hsig  were  separated  into  four  seasons  to  examine  the 
seasonal  influence  on  the  diurnal  variability  (Figure  8).  The  velocity  data  were 
high-pass  filtered  (<  25  hr)  to  enhance  the  diurnal  variability  amongst  the 
seasons..  During  the  spring  months,  the  thermal  gradient  begins  to  increase 
along  the  coast  resulting  in  the  strongest  sea  breeze.  Tsx  has  a  background  of 
0.01  Pa  and  peaks  at  0.11  Pa  at  2300  GMT.  Hsig  increases  0.25  m  to  a  peak  of 
1.55  m  from  a  background  Hsig  of  1.2  m  at  0300  GMT.  In  summer,  thermal 
heating  continues  to  drive  the  sea  breeze  with  Tsx  background  of  0.01  Pa  and 
peaks  at  0.1  Pa  at  2200  GMT.  Due  to  the  eastern  Pacific  Ocean  storm  season 
ending,  the  background  Hs/g  reduces  to  0.9  m  and  diurnally  peaks  at  1.3  m  at 
0300  GMT.  The  summer  season  represents  a  period  when  diurnal  winds  are  the 
most  significant  compared  to  background  waves.  It  also  has  the  largest  effect  on 
generating  sea  beeze  driven  waves,  increasing  Hsig  by  0.4  m.  In  the  fall,  days 
begin  to  get  shorter  resulting  in  a  reduced  coastal  thermal  gradient.  The 
background  Tsx  '\s  0  Pa  and  peaks  at  0.7  Pa  at  2200  GMT.  For  the  first  time,  Tsx 
becomes  slightly  negative.  In  the  fall,  waves  begin  to  be  influenced  by  the  onset 
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of  the  winter  storm  season  with  a  background  Hsig  of  1.15  m  to  a  sea  breeze 
enhanced  peak  of  1.3  m  at  0100  GMT.  In  the  winter  months,  the  thermal 
gradient  continues  to  weaken,  resulting  in  a  background  Tsx  of  -0.01  Pa  with  a 
sea  breeze  peak  of  less  than  0.04  Pa  at  2400GMT.  The  background  Hsig  is  the 
largest  (-1.6  m)  and  no  distinct  peaks  exists,  though  there  is  a  subtle  background 
trend,  which  is  highest  around  0100  GMT  and  slowly  decays  throughout  the  day. 
The  diurnal  forcing  is  the  weakest  in  the  winter. 

Seasonal  surface  and  subaqueous  (below  wave  trough)  cross-shore 
velocities  (Figure  9)  are  examined  as  a  function  of  a  seasonal  diurnal  rsxand  Hsig 
(Figure  8)  The  velocities  in  the  surface  region  are  the  largest  in  the  winter,  owing 
to  larger  waves,  even  though  Tsx  '\s  the  lowest  at  this  time.  The  surface  velocities 
remain  large  in  the  spring,  which  still  have  large  background  waves.  The  summer 
has  the  lowest  surface  velocities  because  this  period  has  the  smallest 
background  waves.  The  surface  velocities  for  all  seasons  vary  diurnally.  The 
diurnal  variability  is  enhanced  in  the  fall,  as  it  a  period  with  relatively  constant 
Hsig,  but  strong  diurnal  Tsx variation. 

For  the  subaqueous  region  (z/h<0.75)  during  sea  breeze  events  (when  Tsx 
>0.05,  2000-2400  GMT),  a  vertically  uniform  offshore  profile  exists.  The  profile 
velocity  magnitude  increases  with  increasing  seasonal  background  wave  height. 
During  the  locally  generated  sea  breeze  waves  (0100-0600  GMT),  there  is 
minimal  enhancement  of  the  velocity  profile,  which  increases  with  increasing 
seasonal  background  wave  height.  During  the  relaxation  period  (0800-1600 
GMT),  the  velocity  profiles  have  vertical  variations  that  vary  seasonally.  Since  the 
Tsx  only  exceeds  the  0.05Pa  threshold  during  the  sea  breeze,  the  velocity  profiles 
vary  primarily  as  a  function  background  waves  over  the  seasons.  The  vertical 
profile  structure  is  consistent  with  subtidal  velocity  structure  (Section  4.1),  Lentz 
et  al.  (2008),  and  Fewings  et  al.  (2008). 
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D.  CROSS-SHORE  TRANSPORT 

The  quantitative  variation  of  surface  drift  as  a  function  of  its  forcing 
parameters  remains  poorly  understood  (Ardhuin  et  al.  2008,  Ardhuin  et  al.  2004, 
Rascle  et  al.  2004).  Ardhuin  et  al.  (2004,  2008)  suggest  that  Stokes  drift  due  to 
wave  forcing  typically  accounts  for  two-thirds  of  the  surface  wind-induced  drift, 
but  these  estimates  are  not  exact  and  require  additional  research.  The  Monterey 
Bay  site  is  dominated  by  a  year-round  sea  breeze  leading  to  a  significant  daily 
wind-induced  drift.  The  theoretical  Stokes  drift  {Qstokes)  is  correlated  with  the 
estimated  surface  transport  {Qsurface.Aocp)  at  the  95%  significance  level  with  a 
slope  (m)  of  0.42  and  a  coefficient  of  determination  (r^  )  of  0.57  implying  that  the 
measured  surface  transport  is  over  estimated  by  -50%.  A  linear  regression  of 
Qstokes  and  Qsurface,ADCP  ,  when  Tsx  was  minimal  (-0.005  Pa  =<  Tsx  <=  0.005  Pa) 
resulted  in  a  m  of  0.89  and  a  of  0.57  suggesting  that  the  difference  in  the 
predicted  and  measured  surface  transports  are  related  to  cross-shore  wind-drift, 
which  is  close  to  the  estimates  by  Ardhuin  et  al.  (2004,  2008)  .  The  difference  in 
the  predicted  and  measured  transports  may  also  be  related  to  the  fact  that  ADCP 
measurements  are  problematic  and/or  the  bin  size  is  not  adequate  to  resolve  the 
transport  in  the  trough-crest  region.  Even  though  potential  biases  exist,  the 
observations  are  used  to  examine  the  diurnal  surface  transport  within  the  trough- 
crest  region. 

Qstokes  is  correlated  with  the  estimated  ADCP  sub-surface  depth  integrated 
transport  {Qsub-surtace.ADCp)  at  the  95%  significance  level  with  a  m  of  0.89  and  an 
of  0.29.  These  results  indicate  a  balance  in  transport  estimates,  but  recognize 
that  this  relationship  is  weak. 

When  the  sea  breeze  develops,  there  is  an  immediate  increase 
(decrease)  in  Qsurtace.Aocp  {Qsub-surface,ADCp)  suggesting  that  the  winds  are 
influencing  the  transport  (Figure  10).  When  the  sea  breeze  induced  Hsig 
increase,,  QsurtaceAocp  and  Qstokes  increases  by  the  same  magnitude  while  Qsub- 
surtace,ADCP  Weakens.  During  the  relaxation  period,  Qsurtace.Aocp  decreases  at  a 
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greater  rate  than  Qstokes,  which  does  not  include  a  theoretical  wind-drift.  Although 
Osurface.ADCP  and  Qstokes  increase  with  roughly  the  same  magnitude  during  the 
diurnal  wave  forcing,  Qsurtace.ADCP  continues  to  be  nearly  double  the  magnitude  of 
Qsfo/ces  during  the  relaxation  period,  which  is  related  to  the  background  wind-drift 
or  ADCP  errors. 

Two  separate  net  transports  were  calculated,  where  the  first  net  transport 
{Qneti)  is  the  sum  of  Qsurtace.Aocp  and  Q  sub-surf  ace,  ADCP  (green  triangles)  resulting  in 
net  onshore  of  5.2  x  10'^  m^/s  and  the  second  net  transport  {Qnet2)  is  the  sum  of 
Qstokes  and  Qsud-surface,/\DCP  (magenta  triangles)  resulting  in  net  offshore  of  5.3  x  10' 
^  m^/s.  During  sea  breeze,  Qnet2  {Qneti)  is  increased  (reduced)  in  the  offshore 
(onshore).  During  relaxation,  Qnef^is  relatively  balanced  in  the  cross-shore,  even 
when  the  sea  breeze  driven  waves  increases,  which  is  consistent  with  the  idea  of 
zero  net  balance  at  depth  (Lentz  et  al.  2008). 

A  slight  rotational  error  in  the  coordinate  frame  of  reference  may  be 
influencing  Qneti  and  Qnet2  ■  As  mentioned  earlier,  the  coordinate  frame 
orientation  was  based  on  a  PCA  of  the  Uda  between  the  normalized  depths  of 
0.25  to  0.75,  which  resulted  in  a  clockwise  rotation  of  the  coordinate  system  by 
10  degrees.  A  subsequent  PCA  of  the  surface  flow  resulted  in  a  clockwise 
rotation  of  16  degrees.  This,  along  with  the  site’s  complex  bathymetry  due  to  its 
proximity  to  the  Monterey  Bay  canyon  suggest  there  may  be  a  slight  discrepancy 
in  the  rotation  coordinate  frame  orientation.  Simple  adjustments  to  the  rotation 
angle  based  on  net  transports  resulted  in  a  required  clockwise  rotation  of  17 
degrees  to  obtain  a  Qneti  of  zero  and  a  required  clockwise  rotation  of  2  degrees 
yielding  a  Qnet2  near  zero.  The  small  changes  in  coordinate  frame  rotation 
required  to  balance  the  net  transports  signify  the  importance  of  a  precise 
coordinate  frame  orientation. 
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V.  SUMMARY  AND  CONCLUSIONS 


Two  years  of  measured  cross-shore  velocity  profiles,  which  varied  as  a 
function  of  wind  and  wave  forcing  were  examined  on  the  inner  shelf  in  Monterey 
Bay  with  particular  emphasis  on  the  diurnal  variability  associated  with  sea 
breezes.  In  the  absence  of  winds,  subtidal  (>25hrs)  cross-shore  velocity  profiles 
consist  of  a  wave-driven  undertow  profile  as  described  by  Lentz  et  al.  (2008)  and 
were  modified  in  the  presence  of  wind  (Fewings  et  al.  2008).  The  diurnal  (<25hrs) 
variability  of  winds,  waves,  and  ocean  currents  represent  approximately  50%  of 
the  total  signal.  On  the  diurnal  (7-25hrs)  temporal  scale,  the  winds  and  waves 
vary  diurnally  with  a  repetitive  nature.  Onshore  winds  increasing  in  the  afternoon 
for  a  short  duration,  referred  to  a  sea  breeze.  During  the  subsequent  land  breeze 
(relaxation),  the  winds  are  reversed  offshore,  but  because  of  the  mean 
background  onshore  wind,  a  true  offshore  wind  stress  rarely  occurs  in  Monterey, 
CA.  The  diurnal  sea  breeze  increases  the  local  waves  with  a  maximum  occurring 
a  few  hours  after  the  maximum  in  the  winds.  When  the  sea  breeze  occurs,  the 
wind  stress  exceeds  0.05  Pa,  which  is  an  arbitrary  threshold  for  when  the  velocity 
profile  is  influenced  by  winds.  The  background  velocity  profile  in  Monterey  is  a 
wave-driven  undertow  profile  that  is  only  modified  during  sea  breeze  events. 
Results  suggest  wind  and  wave  forcing  have  similar  effects  on  the  surface  and 
subaqueous  current  structure  on  both  the  subtidal  and  diurnal  scales.  The 
velocity  profiles  vary  seasonal  owing  to  storms  and  difference  in  thermal 
gradients.  However,  a  similar  diurnal  velocity  profile  trend  exists  with  only  subtle 
seasonal  modifications  associated  with  changing  wind  and  wave  magnitudes. 

The  background  Eulerian  wave-driven  undertow  is  equal  and  opposite  to 
the  Lagrangian  Stokes  drift,  the  net  transport  at  depth  is  zero  resulting  in  no 
cross-shore  transport.  When  the  sea  breeze  occurs,  Tsx  >  0.05  Pa,  the  undertow 
profile  is  modified  and  the  net  transport  at  the  surface  is  shoreward  and  seaward 
near  the  sea  bed,  referred  to  as  a  two  layer  circulation.  During  the  relaxation 

period  (land  breeze),  the  Tsx  is  not  large  enough  to  modify  the  background 
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undertow  profile  or  able  to  reverse  the  previous  sea  breeze  transport.  Therefore, 
material  is  hypothesized  to  incrementally  move  onshore  near  the  surface  and 
offshore  near  the  sea  bed  only  during  sea  breeze  events,  which  are  periodic  at 
the  center  of  Monterey  Bay. 
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APPENDIX  A:  FOURIER  ANALYSIS  ON  A  SIGNAL  COMPOSED 
OF  TWO  SINUSOIDS  WITH  THE  SAME  FREQUENCY  AND 

DIFFERENT  PHASE 


If  a  signal,  g(t),  is  composed  of  two  sinusoids  with  the  same  frequency  and 
a  different  phase,  defined  as 


^(t)  =  ^co5(a>t)  +  Bcos(a>t+  a) 


(a1) 


where  A  and  B  are  the  amplitude  of  the  two  signals,  w  is  the  radian  wave 
frequency,  and  a  is  the  phase  difference.  g(t)  can  be  re-written  using  a 
trigonometric  identity, 

g(t)  =  >lco5(&>t)  -t-  Bcos  a)t  cosa  —  Bsin  cot  sina  ^^2^ 

The  Fourier  in-phase  amplitude,  an,  is  computed  by. 


a„  =  ^  Acos(cot)coscot  +  ^  /J  Bcos(^cot)coscot  cos  a  —  ^  B sina  sinoot  co scot 

(a3) 


2  rT 


2  rT 


resulting  in 


=  A-\-  Bcosa 


(a4) 


The  Fourier  out-of-phase  component,  bn,  is  computed  by 

=  -  /^  Acos(^co€)sin(x>t  +  “  /(,  Bcos(^cot)  sincot  cos  a  ~  ~  Bsina  sincot  sincot 
(a5) 

resulting  in 

i>„  =  Bjimi 

The  Fourier  coefficients  are  composed  of  amplitudes  from  both  sinusoids 
resulting  in  two  equations  with  three  unknowns,  so  without  knowing  the  exact 
phase  difference,  the  two  signals  cannot  be  separated. 
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APPENDIX  B.  FIGURES  AND  TABLES 


Figure  1 .  Bathymetry  view  of  Monterey  Bay,  CA.  Red  star  indicates  the  location 
of  the  ADCP  in  13m  water  depth  and  the  blue  star  indicates  the  location  of 
the  wind  anenometer.  Green  star  indicates  site  location  for  Woodson  et 
al.  (2007). 
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Figure  2.  Principal  Component  Analysis  of  depth-averaged  velocities,  indicating 
a  10°  clockwise  rotation.  The  multiple  color  dots  are  a  scatter  plot  of  the 
twenty-minute  mean  velocities  at  normalized  depths  of  z/h  =  0.25  through 
z/h  =  0.75  at  0.05  increments.  The  blue  lines  indicate  the  major  and 
minor  axes. 
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Cross-shore  Wind  Stress 


Figure  4.  a.)  Twenty-minute  mean  of  Ts  for  2007  b.)  Twenty-minute  mean  of  Ts  for 
July  2007  c.)  Twenty-minute  mean  of  Hsig  for  2007  d.)  Twenty-minute 
mean  of  Hsig  for  July  2007 
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Figure  5.  a.)  Spectral  analysis  of  the  cross  (along)  -  shore  winds  in  blue  (red)  . 
b.)  Spectral  analysis  of  significant  wave  height,  c.):  Spectral  analysis  of 
the  cross-shore  vertical  current  profile,  d.)  Spectral  analysis  of  the 
alongshore  vertical  current  profile.  Calculated  on  a  2  year  record  with  20 
minute  averaged  signal.  Spectra  has  102  DOF  and  A/  =  0.003  cpd. 
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Figure  6.  Hourly  means  of  cross-shore  velocity  in  blue  circles  as  a  function  of 
normalized  ADCP  bin  depth,  a)  when  waves  are  small  (<  0.7m)  and  cross¬ 
shore  wind  stress  is  small  (<  0.05  Pa),  b)  when  waves  are  small  (<  0.7m) 
and  cross-shore  wind  stress  is  large(>  0.1  Pa),  c)  when  waves  are  large  (> 
1.7m)  and  cross-shore  wind  stress  is  small  (<  0.05  Pa),  and  d)  when 
waves  are  large  (>  1.7m)  and  cross-shore  wind  stress  is  large  (>0.1Pa). 
The  numbers  in  the  bottom  right  corner  of  each  plot  represent  the  number 
of  occurrences  of  each. 
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Figure  7.  a.)  Hourly  mean  of  cross-shore  wind  stress  in  Pa.  b.)  Hourly  mean  of 
significant  wave  height,  Hsig  ,  in  m.  c.)  Vertical  profile  of  cross-shore 
current  as  a  function  of  depth. 
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Wind  Stress  (Ra) 


Figure  8.  Seasonal  Variability:  Top  panel)  Spring  (Mar/Apr/May),  Top  Left  panel: 
Hourly  mean  of  cross-shore  wind  stress  in  Pa.  Spring  (Mar-May)  in  red 
circles,  Summer  (Jun-Aug)  in  magenta  stars,  Fall  (Sep-Nov)  in  blue  balls, 
and  Winter  (Dec-Feb)  green  triangles.  Bottom  panel:  Hourly  mean  of 
Significant  wave  height  Hsig  in  m.  Spring  (Mar-May)  in  red  circles. 
Summer  (Jun-Aug)  in  magenta  stars.  Fall  (Sep-Nov)  in  blue  balls,  and 
Winter  (Dec-Feb)  green  triangles. 
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Figure  9.  a.)  Vertical  profile  of  cross-shore  current  as  a  function  of  depth  for 
Spring  (Mar-May),  b)  same  as  for  a.  but  for  Summer  (Jun-Aug).  c)  for  Fall 
(Sep-Nov).  d)  for  Winter  (Dec-Feb).  The  current  data  were  high-pass 
filtered  (cutoff  of  25hr). 
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Figure  10.  Mass  transport;  Qsurface.ADCp  (blue  stars),  Qstokes  (black  squares)  and 
Qsub-sutface,ADCp  (red  circles).  The  first  net  transport  (Qneti)  is  the  sum  of 
Qsurface.ADCp  and  Qsub-surface,ADCp  (green  triangles)  and  is  onshore  at  5.2  x  10' 
^  m^/s.  The  second  net  transport  (Qnet2)  is  the  sum  of  Qstokes  and  Qsub- 
surface  (magenta  triangles)  and  is  offshore  at  5.3  x  10'^  m^/s. 
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Table  1 .  Current  variability  is  the  variance  of  the  high  pass  filtered  (cut  off  of 
25  hr)  velocity  divided  by  the  variance  of  the  total  velocity  at  the 
corresponding  normalized  depths  (z/h).  Tidal  variability  represents  the 
variance  of  predicted  tidal  induced  currents  divided  by  the  variance  of  the 
total  velocity  at  the  corresponding  z/h.  Wind  variability  is  the  variance  of 
the  high  pass  filtered  (cut  off  of  25  hr)  velocity  divided  by  the  variance  of 
the  total  velocity.  Wave  variability  represents  the  variance  of  the  high  pass 
filtered  (cut  off  of  25  hr)  wave  height  divided  by  the  variance  of  the  total 
wave  height. 
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Table  2.  Diurnal  tidal  constituents  (with  a  signal  to  noise  ratio  greater  than 
1.0)  and  their  frequencies  (cph).  Signal  to  noise  ratios  (SNR)  at  the 
normalized  depths  of  z/h  =1.0,  0.75,  0.5  and  0.25  are  calculated  for  each 
diurnal  tidal  constituent  using  the  T_TIDE  program. 
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